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Abstract — In this paper the mean time to get the recoupment of organism due to stress related disease or dysfunction is 
obtained through shock models. Two types of non-cumulative damage shock models are considered here. Based on the 
distribution of damage, caused by a shock affecting a biological organ system, the intervals with small, intermediate and 
large damage are introduced. The initial homogeneous Poisson process is split into three homogeneous Poisson processes 
and studied independently. In this study we also consider the administration of Low salt & High salt diet causes a shock 
effects in terms of aldosterone in the human system. The study reveals the damage caused to the organ system in the 
sense of blood pressure. The study of data from the patients under the influence of high salt food revealed that the higher 
level of stress can cause permanent damage to biological organism. 

Keywords — Aldosterone, Stress, Poisson Process, Shock Model, Hazard rate, Laplace Transform 

AMS Classification— 60E05, 62N 05. 60Gxx, 62Hxx, 62Pxx 

I. INTRODUCTION 

There are many mathematical models to explain how various factors influence the performance of an engineering 
and / or physical system. As these models are well defined the prediction obtained from these models exactly match with the 
real situations. A biological system is being continuously influenced by various biological, environmental and social factors. 
In the absence of exact information about the contributions of such factors to the biological system, developing a model 
which can predict exactly is difficult. The biological activities such as milk production, reproduction etc are greatly 
influenced by factors like heat, cold, body weight, food and other environmental stresses. As it is essential to identify the key 
parameters which affect the production of milk and reproduction capability, it is suggested here a mathematical model which 
could be used to assess the effect of key parameters on the biological activities. 

Human are highly adapt in nature. They are able to fend off stresses, of course to certain limit, which poses a threat 
to an individual's physical and psychological well being. Some of these threats are physical in nature, such as extreme 
temperature or malnutrition's or deficit of water. But for most people, perceived threats to their well being often are 
psychological. For example insecurity of job causes mental agony that affect one's well being. Similarly a failed relationship 
between people triggers distress affecting smooth functioning of a biological system. Regardless of the source of the stress, 
the body responds by activating well-defined physiological systems that specialize in helping a person cope with the stress. 

Biological organisms are usually subject to shocks, which are harmful events that occur randomly in time and 
magnitude, and that can cause a failure or death, respectively. Diseases, viruses, heart attack, stress or more generally, 
demand for energy, can be interpreted as shocks for organisms. The stochastic theory of shocks was extensively studied in 
reliability literature, although there are still a lot of open questions from theoretical and practical points of view. 
Traditionally, two basic cases - the cumulative shock model and the extreme shock model - were considered. The former 
means that the system fails when the cumulative shock magnitude enters some critical region. The latter means that the 
system breaks down as soon as the magnitude of an individual shock goes into some given critical region. 

In this paper the mean time to get the recoupment of stress related disease or dysfunction is obtained through shock 
models. Two types of non-cumulative damage shock models are considered here. Based on the distribution of damage, 
caused by a shock affecting a biological organ system, the intervals with small, intermediate and large damage are 
introduced. The initial homogeneous Poisson process is split into three homogeneous Poisson processes and studied 
independently. In this paper we also consider the administration of Low salt & High salt diet causes a shock interms of 
aldosterone in the biological system. The study reveals the damage caused to the organ system in the sense of blood pressure. 
The concept of shock model and cumulative damage process has been used to determine the expected time to the failure of 
the biological organ system under different assumptions. One can refer the shock model approach in Gurland (1955), 
Bartholomew (1971), Esary, Marshall and Proschan (1973), Barlow and Proschan (1975), Bartholomew and Forbes (1979) 
Shanthikumar & Sumita (1983, 1985), Nakagava T, Kijima M (1989), Ross SM(1996), Lakshmi (2001), Lam & Zhang 
(2004) and Lakshmi & Venkatesh (2009). 
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II. ASSUMPTIONS OF THE MODEL 

1) Stresses are events which cause perturbation to the biological system, leading to its deterioration and consequent 
failures or threatening. 

2) The effect of stress response on the biological system is measured by damage process. 

3) Stress response is the source of hormone secretion. The mineralocorticoid aldosterone is secreted from the adrenal 
cortex in response to stress. 

4) Damages are caused by response of stress at each episode, whose inter arrival times are assumed to be i.i.d random 
variables. 

5) If the total damage exceeds a threshold level, which is itself a random variable, the threat occurs and the biological 
system is recognized as stressed. 

6) The random threshold is assumed due to the fact that the coping ability of stresses varies from one individual 
system to another. 

7) The process which generates stress episode, the sequence of damages and threshold are mutually independent. 

III. MATHEMATICAL MODEL 

3.1 Introduction 

Biological organ systems are subject to a sequence of shocks to be modeled by a stochastic point process 

in [0, oo J . Shocks (or Stresses) are often harmful and can destroy biological system and therefore probability of survival 

under shock (or stress) is of the main interest in this field of biological theory. Each stress response, independent of the 
previous ones, causes a random damage Z with distribution function W(z). This damage is accumulated and when exceeding 
some given level results in failure of a system [1, 2]. The cumulative damage process was studied in many publications [5, 
11]. Assume for simplicity that stress process is the one of the cause of biological organ system failures [7, 12]. 

In the simplest non-cumulative setting each shock independent of the previous ones and the state of a system leads 
to its failure with a fixed probability p. In this case the probability of a system's survival in [0, t) can be easily obtained for 
the Poisson process of shocks. 

It is useful noting that in the alternative interpretation of the non-cumulative model a variable Z is understood not 
like a direct damage to a system but as a random level (shock) of a stress. The probability of a biological organ failure in this 
model depends on the realized value of Z. 

3.2 Splitting the shock process 

A non - homogeneous Poisson process with rate n\t) is assumed to be our model for the process of shocks. 
Consider the following n regions (intervals): 

[0,z 1 ),[z 1 ,z 2 ),....[z n _ 1 ,co) (1) 

The probability that the damage does not exceed level Z, is Wi Z t ) , < I ' < J ' ; i, j <tl , d n = QO, d = is denoted 
by 

P u =w( Zj )-W(z t ) 

P j =w(z j )-W(0)=w(z j ) (2) 

The first important step is to derive the probability P:\t) that all shocks which have occurred in [0, t) had 
resulted in the damage not exceeding C i . This probability can be defined as 

P i {t)=Qy^\-^{i-p i )h{x)dx\ (3) 

Consider a terminating orderly point process when each event with probability 1 — p i it) where t is the time since 
the starts of the process terminates it. Let t denote the random time to termination. The conditional hazard (CH) rate 
A [t, H t ) can be defined as 

A c (t,H t )dt=Pr\T p e[t,t + dt)/H t ,T p (H t )>t\ 

=(l- Pi (t)M?,H t ))dt (4) 
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Where H t and A\t,H t ) denote the history in 10,?) and the complete intensity function of the point process 
respectively. [14] 

The specific case of the Non Homogeneous Hazard process: 

A z (t,H t )=h(x) 

Z z (t,H t )=(l-Pi(t)Wt) (5) 

The probability that all shocks had resulted in damage in the interval \z i , Z ■. ) can be defined as 

P. . (*)=exp - }(l - p ( j )h(x)dx (6) 

Consider the most important specific case, in practice of the three regions [0 , Z s ), [Z s , Z, ) , yZ, , °0 1 where 
"s" stands for "small" and "1" stands for large. In accordance with definition (3), denote probabilities of the damage in the 
corresponding interval as /? s , p s , , p, . Then 

P s (t)=exp l-\(l-p s )h{x)dx\ (7) 



p s j(t)=ew\- ft- Ps,Mx)dx\ 



(8) 



P / (?)=expj-J(l- Pl )h(x)dx\ (9) 

Here P s (t)is the probability that all shocks from the Poisson process in [0, t) result in a "small" damage 
in[0,Zjj; P s i(t) denote the probability that all shocks in an intermediate damage in|z c , Z, ); Pi\t) denote the 
probability that all shocks in a "large" damage in the interval \z t , °° J . 

The strongest non - cumulative criterion of failure in [0, t) is when at least one shock had resulted in a damage 
exceeding Z s • This means that all shocks with damage from the first region are allowed. Then the corresponding survival 

function is defined by P s \tj in Eq. (7). In this stage the stress level is mild. 

Let the shock process be the homogeneous Poisson process with rate h. Then it can be split into three independent 
Poisson processes with rates 

hp s ,hp sl ,hp, (10) 

The operation of splitting [3, 18] is an important tool for the described approach. In accordance with the 
formulated criterion, the failure of our organism will occur in [0,f), if at least one shock from the process with rate hp l will 

occur or, if more than k shocks from the process with rate hp s , will occur (shocks with small damage are allowed 
without restrictions). These considerations lead to the following formula for the survival function: 



i \ k (hp t)' 

p s.k (/)=exp{- hp, t}exp{- hp sl t]^ — — 



(ii) 

!=0 l 

When there is no intermediate region: Z s — Z i we arrive at the formula 

p s ,o (0 = P s (t)= exp {- hp,t} = exp {- h(l - p s >} (12) 
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It is reasonable to assume in the absence of additional information that each shock with damage in I Z c , Z[ ) can still be 
considered as harmless with some probability. This probability can be defined as 

P a (z)=— ,ze[z,,z,) 



It means that when the damage is close to Z , the probability of accepting it is close to 1 and it is close to for the other 

margin. Integrating over all possible values of c gives 

d, 

Pa=\Pa(z) dW{ Z ) 

Then the probability of a system functioning without a failure in 10, t ) for this case 

P s .a (0= e *P (- h Pi ?}exp {- hp sl (l - p a )t\ (13) 

Indeed, the first multiplier in the right hand side is the probability that not a single shock with damage in the third 
region will occur. The second multiplier defines the probability that all shocks (stresses) with damage in the second region 
are classified as acceptable. 

3.3 Time - dependent criterion of failure 

hi our first interpretation of acceptable shock (stress) the biological system was competent of recovering the 
consequences of each shock or stress with a small damage. In a simplified version it was assumed that no time is needed for 
this operation. On the other hand, it is clear that the model can be more adequate, if the time needed for the system to recover 
after the shock is taken into consideration. It is natural to assume that it is a random variable T with a distribution function 
of R(t) (different values of damage need different time of recovering and this fact is described by R(tJ). Thus, if the shock 
occurs while the system still has not recovered from the previous one, then it is a failure (criterion of failure). [13, 19, 20] As 

previously, we want to derive the probability of a failure -free performance in [0, t): P T (t) .Consider the Poisson process of 
shocks with rate A . The following equation for P T [tj takes place: 



dx (14) 



P T (t)= exp {- A t\l + A t)+ f A exp {- Ax}x \a exp {- Ay }R(y)P* (t-x-y)dy 

o L.0 

Where the first term in the right hand side is the probability that there was not more than one shock in [0, i) and the integrand 
defines the joint probability of the following events: 

> the first shock occurred in I X, X + dx) ; 

> the second shock occurred in \X + y, X + y + dy) ; 

> the time between two shocks y is sufficient for recovering (probability R(y)); 

> the organ system is functioning without failures in I X + y, t) 

By P T [t jin equation (14) we denote the probability of system's functioning without failures in [0,0 given that 

the first shock had occurred at t = 0: Thus, it differs from P T [tjby the initial state of the shock process. This is important 
because, if the first shock is fixed at t = 0; then the next one, occurring earlier thanT , will cause a failure. 

Probability P (t) in its turn can be easily obtained from the following integral equation: 

t 

P* (t)= exp (- A t)+ \ A exp (- A x)r(x) P* (t - x)dx (15) 

o 

Interpretation of equation (14) is simpler than of equation (15): the integrand defines the joint probability of the 
first shock occurrence in \X, X + dx) , this time is sufficient for recovering from the first shock at t = 0, and the organ 
system is functioning without failures in [x,t). 

Equations. (14) and (15) can be easily solved by means of the Laplace transform. Denote the Laplace transforms 
of P T [t ), P (t) by P T (s), P ys) respectively. Applying the Laplace transform to equation (15) 
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{s + A)\[-AR{s + A)} 



Where R\S) denote the Laplace transform of R(t). Performing the Laplace transform of both sides of equation. (14) and 
substituting P T [S ) . ; we eventually obtain 

~,s s\i-AR(s + A)\-A 2 R(s + A)+2A 

p M= A — , , 2 ( — ^r — ^T (16) 

{s + A) 2 [1-AR(s + A)\ 

Now we consider the following specific cases of practical importance. 



Case 1: 

Exponentially distributed T '. R{t)=l — exp (— jUt) . 



Then: R (s + A)= 



JU 



(s + A\s + A + ju) 



and 

~ / \ s + 2A+ u 

s + s{2A + ju) + A 

Performing the inverse Laplace transform 

P T (t)= \ exp (s 1 1) + A 2 exp (* 2 1) (18) 

Where S^ , S 2 are the roots of the denominator in equation (17) 

_-(2A + ju)± V(22 + //) 2 - 4 A 2 

and 



s l +2A + ju A s 2 +2A + ju 



A 1 =-i A 



\(1-R( X )) 



Equation (18) gives an exact solution ioiP T (t). In practical applications it is convenient to use simple approximate 

formulas, showing the role of the key parameters of the problem. 
Consider the following reasonable assumption: 

1 
— »T 

A 

(19) 

Relation (19) means that the mean inter-arrival time in the shock process is much larger than the mean time of 
recovery, and this is often the case in practice. In the study of repairable systems the similar case is usually called the fast 

recoup [7, 8]. Using relation (19) for deriving an approximate relation for the probability P T (t) of easel ( T = jU J we 
can obtain P\t)~CXp\—A t) as the second term in the right hand side decreases very sharply with 

t (~exp(-z r ?)). 

The corresponding time dependent error can be easily estimated as 

8 = exp I— A T t 111 — exp I— 2 A T t)\X + oil ) land is usually sufficiently small. The error arising from substitution 

of Aj in Eq. (18) for case 1 is bounded by \At I . The assumption of fast recoup is often considered for highly unfailing 

organ systems ( P (t)is close to 1) [21]. For this case: a ~ 2/t T t . 

Case 2: 

Constant time of recovery: 

, v fO < t < T 

[l t>T 
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General relation will turn to 

~ / \ _ *[l - A exp (- {s + Aft)] - A 2 exp (- (s + Ajr) + 2/t 

(s + A) 2 [l-A exp (- (s + A)f)] (20) 

Here we differentiate above equation with respect to s and setting s=0 we get the expected time to get recovery of stress 
related dysfunction. [12,16] 

d_ ~ / J _ 3 -Acxp(- At\A - At - Aexp(- At)) 
~ds AS l_ - [a{1-c W {-At))Y (21) 

1 

Another and much more effective approximate approach is based on additional assumption (19): — >> T Indeed, 

A 

applying equation similar to Eq. (3) for the specific case of this example 

h(x) = A,p i = q = l-exp(- At)~ At 

In the absence of recovery time (but under the assumption that each shock with probability I — p — p results in the 

failure), according to Eq. (3), the probability of failure free performance in [0, t J is exactly exp (— Aqt ) : 
In the presence of a recovery time the corresponding probability can be written as 

P T (t) = E[exp(-Aq(t-T s ))] (23) 

where T s denote the total random time the system is in the state of recovery in 10, t) . 

Assume that the system is highly unfailing (the probability of failure-free performance is close to 1) and relation (19) holds. 

Using Taylor's expansion in a series in equation (23) for t » E\T s j 

S = exp (- Aqt) - £jexp (- Aqit - t s ))] 

\ | h 4- 2 ) ~ 



A qft 



It 



(1 + (1)) 

(24) 
Relation (24) is derived for an arbitrary R(t). WhenRu) is a step function of Case 2, the term in square 

brackets can be simplified to I + Aqt I Z .It can be easily seen that Eq. (20) follows from equation (24) for the exponential 
case. It also follows from equation (23) that the additional to relation (19) assumption for the definition of the fast recoup is: 



Aq E\t 2 ) 



<< 1 .This means that the second central moment of T should be bounded. But this is the case for the usual 

2? 

lifetime distribution functions to be used for modeling the recovery time. 

IV. APPLICATION 

Effect of salt stress on biological system 

In this study we consider the administration of Low salt (LS) diet & High salt (HS) diet causes a shock in terms of 
aldosterone in the biological system. The study reveals the damage caused to the organ system in the sense of blood pressure. 

Methods 

A group of patients with unilateral aldosterone -producing adenoma (APA) were feeding a LS diet for one week 
followed by HS diet for another week. At the end of each diet period, measurements of blood pressure (BP) at regular 
interval, daily aldosterone measurement and a metaclopramide (MCP) test were performed[21]. 

It was observed that the blood pressure level both systolic and diastolic were remarkably high during the day time 
and night time due to the administration of LS diet. Also it was observed that HS diet caused a steep increase in BP both 
systolic and diastolic in day time and night time. From the observations it was revealed that there was a vast difference 
between systolic blood pressure (SBP) level in day time due to LS diet and diastolic blood pressure (DBP) level in day time 
due to HS diet; but there was no significance in DBP level due to LS & HS diet. The conclusion is the HS diet causes very 
high level of SBP (>140 mm Hg) in a normal person due to the intake of salt diet. 

Results of MCP Test 

Patients diagnosed with unilateral APA were enrolled for MCP test. The patients had abnormally high ratio of 
plasma aldosterone concentration (PAC) to plasma renin activity (PRA), ranging from 52 to 468 with a mean of 189+ 31. 
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After administrating an intravenous injection of 10 mg MCP, the patients PRA was further suppressed by HS. Initially the 
change in PAC was marginal but after few minutes their PAC increased significantly. [6, 22] 
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Fig. 1: Left Panel: Changes in PAC after a 10-mg intravenous administration of metoclopramide on 
low-salt and high-salt. Right Panel: Increment of PAC from basal on low-salt diet and high-salt diet 



V. RESULTS 

The mean increment of PAC at basal level is considered as shock response which is exactly fitted the Poisson 
Process whose rate is h( X ) = A and T is mean time of recovery for the corresponding shock response which is distributed 
exponentially and the realized damage level is measured as in the sense of blood pressure. 



In our application let us take the parameter values A = 1 .45 & T = 0.033 and A = 2.41 86 & T = 0.033 
we get E(T ) = 2.221805 & E(T) = 2.600104 from low salt diet and high salt diet in patients with aldosterone 
producing adenoma respectively. 
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Fig.2: Left Panel: Mean time changes due to aldosterone concentration by using Poisson shock 
model. Right Panel: Survival probabilities of patients on salt manipulation 



VI. CONCLUSION 

Our mathematical model reveals that the negligible amount of stress effects does not affect the organ system 
immediately but if it is too long or it is too often i.e, it survives more time to stand in its positions strongly which induces 
the organ system to failure. In our result, the rate of stress effect between 50 and 120, it gives the small level of damage to 
the system which is mild and also the rate of stress effect between 120 and 130, it causes moderate level of damage to the 
organ system, but the rate of stress effect exceeding above 130, the system fails under this circumstances. In summary, in our 
application part, we defined two subtypes of patients with APA according to their responses to MCP during salt 
manipulation. On HS, the so-called nonsuppressible patients had higher level blood pressure. 

Although different responses of aldosterone secretion to MCP in various hypertensive diseases have been reported, 
the significance of individual variations has not studied. In the present study, we defined two subgroups of patients APA 
according to their different responses of aldosterone secretion to MCP on changes of salt intake. Patients with a larger 
difference in the increment areas of PAC between LS and HS are thought to be more suppressible by Dopamine. Therefore, 
we speculate that angiotensin responsive APA may have a greater response to MCP. Although the absolute PAC after MCP 
was not different between LS and HS, the PAC increment was greater on HS. This indicates that changing salt intake does 
not alter autonomous aldosterone secretion from APA. 
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Finally, we discuss human systems under the influence of two different types of salt stress (LS & HS). Low salt 
stress does not alter PAC then the expected damage level is less than the given threshold level and expected recovery time is 
also negligible, but in the case, high salt induces the PAC level too high than the expected level of damage and the recovery 
time is also high. This paper may be extended to patients of hypokalemia due to hypertension and PRA suppression in all by 
increasing number of subjects and using other medical variables. 

Many mathematical models have been developed for the simulation of various engineering systems such as 
mechanical, electrical and etc. But little progress has been made in constructing mathematical models to simulate various 
biological organ systems. Here we have developed a mathematical model to study the effects of stress in aldosterone. Similar 
approach can be adapted to analyze the effect of hypertension and PRA in patients affected by hypokalemia. Further 
attempts can be made to construct relationships between stresses and psychological effects of biological systems to estimate 
their behavior for varied stress responses. 
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